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Many-body efTects are at the very heart of diverse phenomena found in condensed-matter physics. 
One striking example is the Mott insulator phase where conductivity is suppressed as a result 
of a strong repulsive interaction. Advances in cold atom physics have led to the realization of 
the Mott insulating phases of atoms in an optical lattice, mimicking the corresponding condensed 
matter systems. Here, we explore an exotic strongly-correlated system of Interacting Dual Mott 
Insulators of bosons and fermions. We reveal that an inter-species interaction between bosons and 
fermions drastically modifies each Mott insulator, causing effects that include melting, generation of 
composite particles, an anti-correlated phase, and complete phase-separation. Comparisons between 
the experimental results and numerical simulations indicate intrinsic adiabatic heating and cooling 
for the attractively and repulsively interacting dual Mott Insulators, respectively. 



The Mott insulator transition is a dramatic manifesta- 
tion of quantum many-body effects. In condensed mat- 
ter systems, strong repulsive interactions between elec- 
trons lead to the insulating behavior of materials, when 
the commensurate condition for the density is satisfied. 
Several interesting phenomena also have been discov- 
ered around the Mott transition point, including diverse 
phases showing anomalous fluctuations and orderings in 
the spin, charge, and orbital degrees of freedom. These 
are thought to be deeply connected with the mechanism 
of high-temperature superconductivity in cuprates and 
colossal magnetoresistance in manganites[l|. 

Recently, Mott insulating phases have been studied by 
ultracold atoms confined within an optical lattice 0, Q 
where three mutually perpendicular standing laser waves 
create a periodic potential. Making use of the high de- 
gree of control possible with atomic systems, both the 
bosonic superfluid-Mott insulator transition 1443] and 
the fermionic metal- Mott insulator transitions [8|, |9[ have 
been experimentally demonstrated. The detailed behav- 
ior of the Mott insulating phase has been well studied by 
various spectroscopic techniques SEMI- 

Given the important role of the interaction effects in 
these systems, one could expect drastic changes when two 
different Mott insulators with strong interspecies interac- 
tion are combined together in the same lattice. Such a 
system of strongly interacting dual Mott insulator shown 
in Fig. 1 offers a unique possibility for exploring rich 
phases that are not accessible with a single Mott insu- 
lator system. In particular, in addition to interspecies 
interactions, the relative fillings of two species bring a 
different kind of diversity to this system, which provides 
us with a novel paradigm such as filling induced p/iasespy] . 

Theoretical work in this area has been rich. Various 
novel phases have been theoretically discussed for mix- 
tures of bosons and polarized fermions 1^ Hi , Bose-Bose 



mixtures 17. 18 1, and Fermi- Fermi mixtures [l^- Pio- 
neering experimental work has studied the Bose-Fermi 
mixture of rubidium (^^Rb) and spin polarized potas- 
sium (^°K), where impurity fermions change the bosonic 
superfluid-Mott insulator transition 1^] . However, 
until now, the study of interaction and filling induced 
phases in the dual Mott insulating regime was not ex- 
plored. This is also the case for Bose-Bose mixtures 23- 
25 1 and Fermi- Fermi mixture 2^ 27 1 



In this article, we investigate strongly-correlated sys- 
tems of attractively and repulsively interacting dual Mott 
insulators of bosonic and fermionic atoms. To character- 
ize the dual Mott insulators, we measure double occu- 
pancy of bosonic and fermionic atoms and pair occupancy 
of boson and fermion for varying numbers of atoms. From 
these measurements in combination with numerical cal- 
culation, we reveal how the boson-fermion interspecies 
interaction and the number of atoms drastically modify 
the basic behavior of each Mott insulator, both bosonic 
and fermionic. With repulsive interactions, each Mott 
insulator is phase separated or turns into anti-correlated 
mixed Mott insulating phase depending on the number 
of atoms. In contrast, we observe that in the attractively 
interacting dual Mott insulator various composite par- 
ticles of bosons and fermions are formed. Comparison 
between the experimental results and the numerical sim- 
ulation strongly indicates intrinsic adiabatic heating and 
cooling for the attractively and repulsively interacting 
dual Mott Insulators, respectively. 
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FIG. 1: Schematic illustration of effects of interspecies interaction and filling on dual Mott insulators of bosons 
and fermions. Depending on the sign of the interspecies interaction and filling, the density distributions of each Mott 
insulator of unit occupancy is dramatically altered. The red dashed and blue lines represent the spatial distributions of bosons 
and fermions, respectively. 



THE SYSTEM OF DUAL MOTT INSULATORS 
OF BOSONS AND FERMIONS 



bard Hamiltonian of bosons(i/B) and fermions (-ffp) 
boson-fermion on-site interactions; 



plus 



Our experiment is performed with a quantum degen- 
erate gas of bosonic and fermionic isotopes of ytter- 
bium(Yb) atoms subjected to a three-dimensional opti- 
cal lattice of simple cubic symmetry with a lattice con- 
stant of 266 nm. An attractively interacting Bose-Fermi 
dual Mott insulator is prepared with the bosonic iso- 
tope of ^^°Yb mixed with the fermionic isotope of -'^^^Yb. 
Separately, a repulsively interacting dual Mott insula- 
tor is prepared with the bosonic isotope of -'^^^Yb mixed 
with ^^■^Yb. In both types of dual Mott insulators, the 
fermionic isotope of ^^"^Yb is prepared in an almost equal 
population of six nuclear spin components of / = 5/2. 
The relevant s-wave scattering lengths of Yb bosons and 
fermions used in this work are as follows: ai7o=3.38 
nm, ai74=5.55 nm and ai73 = 10.55 nm as well as the 
interspecies s-wave scattering lengths, 0170-173=— 4.30 
nm and 0174-173=7. 34 nm. Here, it is noted that 0173, 
ai7o-i73 and 0174-173 are independent of nuclear-spin 
components of ^'''^Yb|28l|. 

The physics of a Mott insulator is well described by 
a single-band Hubbard model. It considers particle- 
tunneling between adjacent lattice sites and on-site in- 
teraction. In optical lattice experiments, the underlying 
harmonic confinement gives rise to an additional spatially 
dependent energy offset. This has important influence on 
the spatial structure of the Mott insulators f29']. 

As a result, our system of dual Mott insulators of 
bosons and fermions can be described by the Hub- 



HbF —Hb + Hf + UbF flB,i'nF,i,cT, 



Ubb 



+ ^ ^iflB.i, 



nB,i(nB,i ~ 1) 



(1) 



(2) 



i,(7 



where denotes summation over nearest-neighbors 

and a G {1,2, ...,6} denotes the spin component of 

fermionic atoms. 6| (//o-) ^^'^ ifi,(y) are creation 
and annihilation operators for bosons (fermions with a 
spin). nB,i = blbi (np.i.o- = flafi,cr) are the number 
operators for bosons (single spin component fermions). 
C^BB, C^FF, and Ubf represent on-site interaction ener- 
gies between bosons, fermions with different spin com- 
ponents, and boson-fermion pairs, respectively. Here 
Ubb < \Ubf\ < Uff is satisfied in both systems. J is 
the tunneling matrix element, and is the energy offset 
due to an external confinement potential. 



3 



DOUBLE AND BOSON-FERMION PAIR 
OCCUPANCY MEASUREMENTS 

In order to investigate the dual Mott insulators, we 
measure site occupancies, such as double occupancy of ei- 
ther bosons or fermions and pair occupancy of boson and 
fermion, which are schematically depicted in Fig. 2(a). 
The suppression of double occupancy is a key feature of 
the incompressibility of a Mott insulator Q and the pair 
occupancy of boson and fermion is a direct measure of 
their overlap. 

To measure the number of double occupancy, we use a 
one-color photoassociation (PA) technique. For the site 
with two atoms, the irradiation of the PA laser beam 
results in the formation of excited state molecule which 
escapes from the trap, which can be seen through atom 
loss. We apply a strong enough PA laser beam so that 
all the doubly occupied sites are depleted with high ef- 
ficiency. As occupancy of more than two bosons or two 
fermions is negligible in the current experimental param- 
eters, the atom loss corresponds to the absolute number 
of double occupancyjs^. Figure 2(b) shows the atom 
loss due to a strong PA laser beam where the number of 
remaining atoms is obtained from absorption images. 

An illustrating example of PA spectra is shown in Fig. 
2(c) with bosonic homonuclear resonances in an opti- 
cal lattice. In the case of a pure bosonic cloud with a 
Mott plateau of one atom per site, there is essentially 
no atom loss, as expected. When we add fermions to the 
bosonic Mott insulator, however, we observe a significant 
atom loss, meaning sizable double occupancy of bosons. 
This suggests that a fermion admixture induces a dras- 
tic change in the bosonic Mott insulating phase. This 
novel phenomenon of fermion-induced bosonic double oc- 
cupancy is discussed in detail in the following sections. 

The pair occupancy of boson and fermion is measured 
in the same manner by utilizing heteronuclear PA res- 
onances between bosonic and fermionic isotopes. The 
detailed information on the resonances is given in Meth- 
ods. 
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REPULSIVELY INTERACTING DUAL MOTT 
INSULATORS 

First, we discuss the case of repulsively interacting dual 
Mott insulators. The results of the measurements are 
shown in Fig. 3 (a) , where the bosonic double occupancy 
(Db), the fermionic double occupancy (Dp), and boson- 
fermion pair occupancy (Dbf) are plotted as a function 
of the number of fermionic atoms (N-p). The numerical 
calculations for density distributions for bosons ((nB,i)), 
fermions {{fiF.i}) and the bosons and fermions (('t-b,! + 
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FIG. 2: Photoassociation and occupancy measure- 
ments in an optical lattice, (a) Schematic illustration of 
three types of occupancy in an optical lattice. A bosonic atom 
is represented by red, and fermionic by green. The number of 
atoms circled with a solid line is measured in the respective 
measurement, (b) Decay of the number of bosonic atoms due 
to the irradiation of the PA laser in the bosonic Mott insula- 
tor. Only the doubly occupied sites are depleted by the PA 
laser, (c) PA resonance spectra of ^'^^Yb bosons in an at- 
tractively interacting Bose-Fermi system('^^''Yb-^^^Yb, blue 
circles) and in a pure bosonic cloud that forms Mott plateau 
of unit fiUingC^^^Yb, red squares). The baseline of the spec- 
trum of the pure bosonic cloud is shifted vertically by 1 x 10^ 
for clarity. The solid blue line is a fit with a Lorentzian curve, 
and the solid red line represents no resonance. 



np.i)) are shown in Fig. 3 (c) for various numbers of 
fermionic atoms. The calculated distributions of double 
occupancies for bosons ((Z^B.i)) and fermions ((Dfj)) are 
shown in Fig. 3 (d), and those of boson- fermion pair 
occupancy ((-DBF.i)) in Fig- 3(e). 

From these results we have found that the system is 
characterized by the two distinct phases depending on 
A''^. Here we begin by explaining the calculated results. 
As shown in Fig. 3(c), when Np < IC*, the density 
distributions of bosons ((nB,i)) and fermions ((riF.i)) ap- 
parently overlap with each other, while the total density 
distribution {{fiB.i + nF.i)) clearly satisfies nB,i,-|-nF,i=l, 
meaning the formation of the Mott plateau at the unit 
filhng. In addition, from Fig. 3(d) and (e), bosonic, 
fermionic and boson-fermion pair occupancies take values 
near zero. We thus identify this phase as a novel kind of 
mixed Mott state of bosons and fermions by noting that a 
total set of bosons and fermions satisfies comensurability, 
but species separately does not. 

As the number of fermionic atoms increases, a dif- 
ferent phase appears. Figure 3 (c) clearly shows that 
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FIG. 3: Occupancy measurements and numerical simulation of occupancy distributions. (a)Measured pair occu- 
pancies for the repulsively interacting system; (a) (top) Pair occupancy of bosonic atoms (Db) for A^b = 6 x 10'^ represented by 
red circles, (a) (middle) Pair occupancy of fermionic atoms (Dp) for A'b = 9 x 10'' represented by red circles, and the fermionic 
pair occupancy for pure fermionic atoms (yellow circles) shown for comparison, (a) (bottom) Boson-fermion pair occupancy 
(Dbf/A^b) for iVB = 6 X 10^ (red circles), A^b = 8.5 x 10^ (orange circles), JVb = 13 x 10^ (green circles), A'^b = 15 x 10^ (blue 
circles), and A''b = 20 x 10'^ (purple circles). The data with almost the same mean trap frequency but with slightly different 
geometry (27rx(150, 30, 150) Hz) are represented by blue diamonds (TVb = 5.5 x 10^). Note that the pair occupancy is divided 
by the number of bosonic atoms to show fractional overlap between bosons and fermions. (b) Measured pair occupancies for 
the attractively interacting system; (b)(top) Pair occupancy of bosonic atoms for A'b = 5 x 10^ represented by blue circles, 
(b) (middle) Pair occupancy of fermionic atoms for A'b = 5 x 10^ represented by red circles, (b) (bottom) Boson-fermion pair 
occupancy divided by the number of bosonic atoms (Dbf/Nb) (red circles). The solid lines in (a) and (b) correspond to the 
numerical simulation for (a) rini=35 nK and (b) rini=45 nK where Tini is the initial atom temperature before loading into 
an optical lattice. The numerical results with different initial temperatures by ± 5 nK are also represented by the shaded 
areas along the solid lines. (c)-(e)Numerical simulation for the repulsively interacting system with A^b = 6 x 10^, rini=30 
nK; (c)Column density distributions for bosons ((riB.i)) (red dashed lines), fermions ((nF,i))(blue lines), and their summation 
((nB,i -I- nF,i))(green hues), (d) On-site double occupancy distributions of bosons ((-Db.i)) (red dashed lines) and fermions 
{{DF,i)) (green lines), (e) On-site boson-fermion pair occupancy distribution ((Dbf,!))- (f)-(h)Numerical simulation for the 
attractively interacting system with A'^b = 5 x 10^, rini=40 nK; (f) Column density distributions of bosons ((nB,i)) (red dashed 
lines) and fermions ((nF,i))(blue lines), (g) On-site double occupancy distributions of bosons ((Dbj)) (red dashed lines) and 
fermions ((Df,!)) (green lines), (h) On-site boson-fermion pair occupancy distribution ({Dbf.i))- 
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bosons and fermions become phase-separated for larger 
numbers of fermions. The boson-fermion pair occupancy 
also takes a value near zero, even at larger numbers of 
fermions, as shown in Fig. 3(e). This phase-separation 
can be understood from the relative strength of the on- 
site interaction parameters {Ubb < Uqf < Uff)- Since 
the fermion on-site interaction Uff is the largest, it is 
energetically favorable for fermions to spatially extend 
to avoid double occupancy. Due to the Bose-Fermi re- 
pulsive interaction Ubf, bosons and fermions repel each 
other, which results in a phase-separation. Figure 3(c) 
also shows that bosonic double occupancy grows in the 
center of the trap. 

The experimental results of the occupancy measure- 
ments (Fig. 3(a)) confirm the above-mentioned features 
of numerical simulation, demonstrating that these phases 
are realized in our experiments. First, the bottom panel 
of Fig. 3(a) shows strong suppression of Dbf, lower than 
10 % of A^B, which is consistent with the numerical sim- 
ulation (Fig. 3(e)). Df stays small, as shown in the 
middle panel, which is consistent with (-Dp.i) close to 
zero(see Fig. 3 (d)). The top panel shows the suppres- 
sion of Db in the mixed Mott regime at small numbers 
of fermions, which is consistent with the numerical cal- 
culation of {DB,i) in Fig. 3 (d). At larger numbers of 
fermions, we observe a large Db in Fig. 3 (a). We can 
explain this phenomenon by the strong compression ex- 
erted on bosons from the surrounding fermion shell, due 
to the strong interspecies repulsion. Numerical calcula- 
tions using the temperature as the single fitting parame- 
ter quantitatively explain all the experimental data. 
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FIG. 4: Various types of composite particles in the 
attractively interacting system. The on-site interaction 
energy of the composite particles depends on the composition 
of particles. Type A, B, C, and D composite particle have 
energies of ?7bb -I-2[7bf = 1.19[/bf, Ubf, C/bb + ?7ff -|-4f/BF = 
0.67C/bf, and Uff + Ubf ~ —0.63Ubf, respectively. 



ATTRACTIVELY INTERACTING DUAL MOTT 
INSULATORS 

A remarkable difference between the repulsively and 
the attractively interacting system is observed. This 
highlights the crucial role of interspecies interaction in 
strongly interacting dual Mott insulators of bosons and 
fermions. 

The experimental results are summarized in Fig. 3 (b). 
The clear departure from the behavior of the repulsively 
interacting system is seen in the boson-fermion pair occu- 
pancy measurement, where we observe the significantly 
large pair occupancy. The behavior of bosonic double oc- 
cupancy, shown in the top panel, is also quite different. 
The measured Dg shows a characteristic nonmonotonic 
dependence on Nf values and three phases can be iden- 
tified. 

In phase I, where Nf is smaller than 2.5 x lO'', Db 
rapidly increases as Np increases. In the regime of inter- 
mediate number of fermionic atoms (phase II), Db grad- 
ually decreases, while Dbf gradually increases. In the 
regime of large number of fermionic atoms (phase III), 
Db increases again and Df also increases to a significant 
level. 

Each phase is characterized by a formation of differ- 
ent types of on-site composite particles 16| near the trap 
center. The transition from one phase to another can 
be understood qualitatively by considering the relative 
strength of the various types of on-site interactions and 
the associated composite particles, as shown in Fig. 4. 

In phase I, a composite particle of two bosons and one 
fermion (designated type A ("BBF")), which has the low- 
est on-site interaction energy, is predominantly formed 
until the number of fermions reaches the half of the num- 
ber of bosons. This phase is thus characterized as a rapid 
increase of Db and Dbf , as is observed in our experiment 
(Fig. 3 (b), top and bottom panels). 

In phase II where the number of fermionic atoms is 
intermediate {Nf > 2.5 x 10'^), singly occupied fermions 
should appear because only a limited number of fermions 
can form type A particles ("BBF"). Note that it is more 
energetically favorable to form two composite particles of 
one boson and one fermion (designated type B ("BF")) 
than to form one singly occupied fermion and a type 
A particle ("BBF"). Therefore, upon an increase of the 
number of fermions, already formed type A particles 
( "BBF" ) are dissolved into two type B particles ( "BF" ) . 
Namely, "BBF"-f "F" "BF"-f "BF". This results in a 
decrease of bosonic double occupancy and increase in 
boson-fermion pair occupancy, which is again consistent 
with our observations (Fig. 3 (b), top and bottom pan- 
els). 

In phase III, where the number of fermionic atoms is 
large, Df begins to grow in the center due to the har- 
monic confinement. It is important to note that a com- 
posite particle of two bosons and two fermions (desig- 
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nated type C ("BBFF")) is energetically favorable than 
the case of one singly occupied boson and a composite 
particle of one boson and two fermions (designated type 
D ("BFF")), as well as the case of two bosons, singly or 
doubly occupied, and one doubly occupied fermions. As 
a result, the stair- wise density distribution of fermionic 
atoms attracts bosonic double occupancy to the center 
and type C particles ("BBFF")) are formed. Namely, 
"BFF" + "B" "BBFF". This is again consistent with 
the observed increase of fermionic and also bosonic dou- 
ble occupancy shown in phase III of Fig. 3(b). 

Numerical calculations quantitatively verify the above 
described models. In fact, the density distributions 
of bosons ((nB,i)) and fermions ((nF,i)) (Fig. 3 (f)), 
and double occupancies of bosons ((-DB,i)) and fermions 
((^B,i)) (Fig. 3 (g)), and the boson-fermion pair occu- 
pancy (DBF.i) (Fig. 3 (h)) well reproduce the overall 
behaviors discussed above. Note that the adiabatic heat- 
ing, which we discuss in the next section, takes place and 
the boundary between each phase is rather smoothed. 



THERMODYNAMICS 

As is discussed in the previous sections, the numerical 
calculations show reasonable agreement with the experi- 
mental data within a small uncertainty of the initial tem- 
peratures Tini = 40 ± 10 nK. Note that the loading into 
the lattice is assumed to be adiabatic and thus isentropic, 
and that the drastic change occurs in the atom tempera- 
ture during this adiabatic loading process. Here we care- 
fully evaluate the final temperatures after lattice loading, 
^lattice) through the comparison between the experimen- 
tal data and the numerical calculations, and reveal the 
adiabatic heating and cooling effects in the attractively 
and repulsively interacting dual Mott insulators, respec- 
tively. The detailed procedures of our numerical calcula- 
tions are described in the Method section. 

In Fig. 5(a) and (b), we re-plot the measured bosonic 
double occupancies along with the theoretical curves for 
several values of Tiattice in the repulsively and attrac- 
tively interacting systems, respectively. We determine 
^lattice realized in the experiments by comparing these 
measured Db values with calculated ones, and the re- 
sults are summarized in Fig. 5(c) and (d). We see that 
^lattice ranges from 5 nK to 40 nK and from 40 nK to 100 
nK for the repulsively interacting and the attractively in- 
teracting cases, respectively. Since Tini is roughly 40 nK 
in our experiments, the results in Fig. 5(c) and (d) show 
that cooling and heating occur in the repulsively and at- 
tractively interacting dual Mott insulators, respectively. 
Although three-body losses and non-ideal conditions in 
real experiments may affect the estimated Tiattico, the 
obtained agreement between the estimated Tiattice (filled 
circles in Fig. 5(c) and (d)) and calculated Tiattice un- 
der the assumption of adiabatic loading with Tini=40±10 



nK (shaded area in Fig. 5(c) and (d)) suggests that the 
attractively (repulsively) interacting system is adiabati- 
cally heated (cooled). 

Here we discuss the mechanism of adiabatic cooling 
and heating. In the repulsively interacting system where 
mixed or phase-separated Mott insulating states appear, 
the spin degrees of freedom in fermionic ^^^Yb can carry 
a large entropy of In 6 per site in the Mott insulator. 
Especially in the mixed Mott states, bosonic atoms can 
also contribute to the entropy as the seventh degree of 
freedom in addition to the fermionic atoms, resulting in 
a further large entropy of In 7. Therefore, under the isen- 
tropic condition, the adiabatic lattice loading towards the 
Mott regime results in the cooling of the other degrees of 



freedom |31l-l34j. Note that the state of phase-separated 
Mott insulators is easily destroyed by thermal fluctua- 
tions, as discussed in Supplementary Information, and 
Db in the phase II rapidly decreases with the increase in 
Tiattice as shown by the theoretical curves in Fig. 5(a). 
Thanks to the cooling effect caused by the large entropy 
carried by ^^^Yb fermions, we can successfully find the 
phase-separated state. 

On the other hand, in the attractively interacting 
system, the system is adiabatically heated over a wide 
range 35]. As already discussed in the previous section, 
a key feature of the attractively interacting system is 
the formation of various types of composite particles, de- 
pending on the number of fermionic atoms. The impor- 
tant consequence of the formation of composite particles 
is the reduction of the entropy due to the reduction in 
the number of particles involved in the system. Thus, 
the adiabaticity of the lattice loading results in the heat- 
ing of the other degrees of freedom and in the melting of 
the Mott plateau. We confirm the correlation between 
Tiattice and the fraction of the composite particles (see 
Supplementary Information). 



OUTLOOK 

We have investigated quantum phases of dual Mott 
insulators of bosons and fermions focusing on the roles 
of interspecies interactions and relative fillings of two 
species. By carefully comparing measured data with the- 
oretical simulations, we have successfully clarified the 
characteristics of several quantum phases and the non- 
trivial thermodynamics during the lattice loading. The 
present work paves the way for exploring diverse quan- 
tum phases around the Mott transitions at lower temper- 
atures such as the alternating Mott insulator |l5l| and the 
novel phases starring composite particles [16|. For fur- 
ther investigations, we have found the important similar- 
ity between our finding and a basic problem in condensed 
matter systems. In some manganates, it is known that 
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FIG. 5: Thermodynamics of repulsively and attractively interacting dual Mott insulators. (a),(b) Bosonic pair 
occupancies for several temperatures after the lattice loading, Tiattice, in (a) repulsively and (b) attractively interacting dual 
Mott insulators. The solid lines represent the calculations, and the circles are the experimental data, (c) Tiattice evaluated 
through comparing the calculations and the measurements shown in Fig. 5(a) (red circles) and calculated Tiattice for Tini=30 
nK(blue dotted line), 40 nK(blue solid line), and 50 nK(blue dash-dotted line), (d) Tiattice evaluated through comparing the 
calculations and the measurements shown in Fig. 5(b) (blue circles) and calculated Tiattice for Tini^SO nK(red dotted line), 40 
nK(red solid line) and 50 nK(red dash-dotted line). 



the metal-insulator transition accompanies the ferromag- 
netism [1] . We observe in the repulsively interacting sys- 
tem that bosonic atoms undergo the crossover between 
the compressible and incompressible states accompanied 
with phase separation. This can be understood as a man- 
ifestation of such a ferromagnetic-type transition, which 
causes the phase separation in the presence of confine- 
ment. This feature provides a new concept of ferromag- 
netism in the strongly correlated many-body systems. 



METHODS 

Preparation of dual Mott insulators of bosons and 
fermions with Yb isotopes 

A quantum degenerate Bose-Fermi mixture of repul- 
sively interacting ^''^Yb-^^'^Yb or attractively interact- 
ing ^^Yb-^^'^Yb is first created in an optical trap pro- 
duced by crossed laser beams with a wavelength of 532 
nm by sympathetic evaporative cooling of bosonic and 



fermionic species |36| . While we prepare the mixtures 
with various numbers of atoms up to 5x10*, the final 
potential depth is always set to the same value. The 
final trapping frequencies of the crossed optical dipole 
trap are 27rx(125, 30, 190) Hz, which are measured by 
a parametric excitation and a dipole oscillation, and the 
temperature is 40 ± 10 nK. A ^™Yb single species Bose- 
Einstein condensate(BEC) is prepared by blasting ^''■^Yb 
with a laser beam resonant to the "'^So ^ '^Pi transition 
after the production of the quantum degenerate mixture, 
and subsequently waiting for thermalization for typically 
100 ms to 250 ms. Single species quantum degenerate 
gases of ^^*Yb BEC or ^^■^Yb Fermi degeneracy are ob- 
tained by a straightforward evaporative cooling without 
resort to the sympathetic cooling. After preparing the 
quantum degenerate mixture in the optical dipole trap, 
we load the atoms into a three dimensional optical lat- 
tice with a lattice constant of 266 nm in a simple cubic 
symmetry. We ramp up the lattice potential smoothly 
to the final value within 200 ms without changing the 
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optical dipole trap potential. We set the final lattice 
depth to 15 Es, for the repulsively interacting system 
and 18 for the attractively interacting system. Here, 
E^ = h? /{2m\^) ~ /i X 4.05 kHz is the recoil energy of 
the lattice laser beam, m is the mass of ^''^Yb and A(= 
532 nm) is the wavelength of the lattice laser. The trap 
frequencies of external harmonic confinements in the op- 
tical lattice are 27rx(145, 81, 206)Hz and 27rx(147, 84, 
209)Hz for 15Er and 18Er, respectively. The laser in- 
tensity of each lattice axis is precisely controlled by the 
feedback of the monitored laser power. The lattice depth 
is calibrated by observing the diffraction patterns due to 
a pulsed lattice beam. 

Here it is noted that one of the features of using Yb iso- 
topes is the negligible difference in the gravitational sags 
which comes from the fact that the mass and trapping 
potentials are almost the same for Yb isotopes. There- 
fore, this clearly focuses the physics on the role of inter- 
species interaction between bosons and fermions in dis- 
cussing the behaviors of our dual Mott insulators. It is 
noted that on-site interaction energies between fermions 
and those between boson and fermion are independent 
of the fermionic spin components, realizing enlarged spin 
symmetry in the system 27|. In this work, we espe- 
cially explore strongly-correlated regime in an optical 
lattice, where interaction energies dominate the tunnel- 
ing energies or thermal fluctuation. Due to the differ- 
ent s-wave scattering lengths, bosonic isotopes of ^^"'Yb 
and ^^°Yb enter the Mott insulating regime of unit fill- 
ing for lattice depths deeper than 11.2i?R and 12.9i?R, 
respectively [33, ISg] . 



Double and pair occupancy measurement 

We use one-color PA resonances associated with the in- 
tercombination transition (^Sq Sq O +^ Pi). The 
PA laser intensity and pulse duration are typically 0.5 '^2 
W/cm^ and 10~30 ms, respectively, which assures of the 
negligible photon scattering loss for the singly occupied 
sites due to the PA light. In addition, to freeze out the 
site occupancy during the PA process, we ramp up the 
lattice depth to 25 E^i in a short time of 2 ms which is at 
the same time long enough to suppress the non-adiabatic 
excitation to higher bands. The PA resonances used for 
the homonuclear cases of ™Yh, ^^^Yb, and ^^'^Yb are lo- 
cated below the corresponding atomic transitions by 618 
MHz, 687 MHz, and 796 MHz, respectively. The PA res- 
onance for the heteronuclear case of -'^''^Yb-^^^Yb {™Yh- 
"3Yb) is located below the ^Sq o ^P^ of ^''^Yb (i™Yb) 
atomic transition by 674 MHz (795 MHz). To suppress 
Zeeman splittings of the PA resonance due to spin de- 
grees of freedom in ^'''^Yb, we reduce the residual mag- 
netic field less than 50 mG by performing high-resolution 
spectroscopy associated with the magnetic sub-levels in 
the metastable •^P2 state of ^''"'Yb. 



Single-band Hubbard model 

We determine the parameters included in the Bose- 
Fermi Hubbard Hamiltonian, i.e., Ub,Uf,Ubf, and so 
on, based on the Bloch wave functions obtained in the ab 
initio manner. From the preliminary calculations con- 
sidering the higher Bloch bands, we have confirmed that 
the contribution of the higher Bloch bands can be neg- 
ligible owing to the low filling of bosonic and fermionic 
atoms in the present experiments. As a result, multior- 
bital effects, broadening effects of the Wannier functions 
[lot , and self-trapping effects [11], hardly emerge in our 
experiments. 

Some experimental parameters such as atomic mass 
and optical dipole trap potential take almost the same 
values for Yb isotopes, which significantly simplifies the 
model. Thus the tunneling matrix elements Jb and Jp 
and the energy offset due to external confinement po- 
tential Ci are set to be equal for bosons and fermions in 
our theoretical analysis. We fix the single-band Hubbard 
parameters with the values relevant to the experiments: 
C/bb ^hx 2.6kHz, C/bf = ^ x 3.5kHz, [/pF = ft- x 5.0kHz 
and J = hx 0.026kHz for the repulsively interacting sys- 
tem, and Ubb = h x 2.0kHz, Ubf — hx —2.5kHz, Uyf — 
h X 6.2kHz and J ~ h x 0.014kHz for the attractively 
interacting system. 



Details of Numerical Simulation 

We quantitatively analyze the large three-dimensional 
optical lattice system at finite temperatures. Here we 
summarize our theoretical procedure. We deal with the 
bosonic Hamiltonian based on the Gutzwiller type ap- 
proximation [40| : 



Hf = -Y. Jf^ + H.c. + ^ E "B.("B. - 1 

i i 



(4) 



where 



Jf = j 



j'£ncighbors 



(5) 



In addition, since the interaction between fermionic 
atoms is stronger than the bosonic interaction, we neglect 
the hopping term of fermionic atoms. We thus obtain 
the effective local Hamiltonian. Then, thermal states of 
bosons and fermions at a certain temperature Tiattico are 
calculated self-consistently such that the expectation val- 
ues {bj) converges, using exact diagonalization method by 
appropriately eliminating the states consisting of a large 
number of bosons. Here, Tiattice is the final temperature 
after loading atoms into the lattice. Such a single-site 
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approximation becomes fairly good when we focus on a 
deep lattice as in this work. 

Currently there are no reliable methods to measure 
TiatticG in experiments. We can instead evaluate this 
Tlattce numerically with the assumption that the lattice- 
loading process is adiabatic and thus isentropic |3ll, |32. 
35l . |4l| . It is easy to calculate the entropy S'lattice of the 



system described by the effective Hamiltonian mentioned 
above by the exact diagonalization method. The initial 
entropy of trapped atomic gases without the lattice, 5'ini, 
is simply given by neglecting the interactions between 
trapped atoms: 



S,,,,/kB=n'Nj.'^ + 3.6NB(^ 
-If \ Jc 



(6) 



where Tini is the initial temperature of the trapped 
atomic gas before lattice ramping up, Tp is the Fermi 
temperature, Tc is the BEC transition temperature, and 
is the Boltzmann constant. 
This approximation of neglecting atomic interactions 
may cause some deviations of our numerical calculations 
from the experimental results. We therefore consider Tini 
as a unique fitting parameter in our simulations, and in- 
troduce small uncertainty of Tini (±10 nK) so as to in- 
clude certain numerical artifacts caused by these theo- 
retical approximations. In addition, it has been pointed 
out that, in the case of mixtures consisting of bosons and 
polarized fermions, the entropy of atomic gases without 
the lattice hardly depends on the interspecies interactions 

It is useful to explain how we analyze the photoas- 
sociation (PA) loss measurements employed in the ex- 
periments. These measurements roughly provide us 
with the information on bosonic, fermionic, and boson- 
fermion pair occupancies at each site. However, in the 
present system consisting of bosons and multi-component 
fermions, complex configurations of atoms are realized at 
finite temperatures. Therefore, the atom loss due to PA 
process is given by the statistical expectation of the fol- 
lowing local operators that count respectively the num- 
ber of pairs corresponding to the bosonic and fermionic 
homonuclear and the heteronuclear PA measurements: 

Z)b,j 2 ^ &00T{nB,i/2)\nB,i,nF,i){nBA,nF,i\, 



(7) 



(8) 



Dbf,!^ ^ min(n_F_i,nB,i)l'^-B,i,n-F,i)(?^B,i,?^_F,i|, 

TLB i,nF i 

(9) 

where riF,i) represents a local Fock state with nB,i 

bosons and nF,i fermions, and nF,i = '^„nF,i.a- Here, 



floor(a-) is the floor function which returns the largest 
integer not greater than its input, and min(a;,y) returns 
the minimum inputs. We straightforwardly calculate the 
sums of expectation values of these local operators over 
the lattice sites, and compare them with the experimen- 
tally measured PA losses. In addition to the above we 
further measure the heteronuclear PA loss after remov- 
ing bosonic double occupancy using homonuclear PA as 
shown in Supplementary Information. The local operator 
of this measurement is given by 

DB-BF^^ = (10) 

min(nF,i, mod 2) \nB,i, nF,i){nB,i,nF,i\, 
where, mod is the modulo function. 



[1 
[2: 

[s: 

[6: 

[7: 
[s: 
[9: 
[10: 

[11 

[12 

[13 
[14 



Electronic address: 'ssugawa®scphys. kyoto-u.ac.jp' 
Imada, M., Fujimori, A. & Tokura, Y. Metal-insulator 
transitions. Rev. Mod. Phys. 70, 1039-1263 (1998). 
Jaksch, D. & ZoUer, P. The cold atom Hubbard toolbox. 
Ann. Phys 315, 52-79 (2005). Special Issue. 
Bloch, I., Dalibard, J. & Zwerger, W. Many-body physics 
with ultracold gases. Rev. Mod. Phys. 80, 885-964 
(2008). 

Fisher, M. P. A., Weichman, P. B., Grinstein, G. & 
Fisher, D. S. Boson localization and the superfluid- 
insulator transition. Phys. Rev. B 40, 546-570 (1989). 
Jaksch, D., Bruder, C., Cirac, J. I., Gardiner, C. W. & 
ZoUer, P. Cold bosonic atoms in optical lattices. Phys. 
Rev. Lett. 81, 3108-3111 (1998). 

Greiner, M., Mandel, O., Esslinger, T., Hansch, T. & 
Bloch, I. Quantum phase transition from a superfluid to 
a Mott insulator in a gas of ultracold atoms. Nature 415, 
39-44 (2002). 

Spielman, I. B., Phillips, W. D. & Porto, J. V. Mott- 
insulator transition in a two-dimensional atomic hose gas. 
Phys. Rev. Lett. 98, 080404 (2007). 

Jordens, R., Strohmaier, N., Giinter, K., Moritz, H. & 
Esslinger, T. A Mott insulator of fermionic atoms in an 
optical lattice. Nature 455, 204-207 (2008). 
Schneider, U. et al. Metallic and insulating phases of 
repulsively interacting fermions in a 3d optical lattice. 
Science 322, 1520 (2008). 

Campbell, G. K. et al. Imaging the Mott insulator 
shells by using atomic clock shifts. Science 313, 649- 
652 (2006). 

Foiling, S., Widera, A., Miiller, T., Gerbier, F. & Bloch, 
I. Formation of spatial shell structure in the superfluid 
to Mott insulator transition. Phys. Rev. Lett. 97, 060403 
(2006). 

Gemelke, N., Zhang, X., Hung, C.-L. & Chin, C. In situ 
observation of incompressible Mott-insulating domains in 
uhracold atomic gases. Nature 460, 995-998 (2009). 
Bakr, W. S. et al. Probing the superfluid-to-Mott insu- 
lator transition at the single-atom level. Sctence 329, 
547-550 (2010). 

Sherson, J. F. et al. Single-atom-resolved fluorescence 
imaging of an atomic Mott insulator. Nature 467, 68-72 



10 



(2010). 

[15] Titvinidze, I., Snoek, M. & Hofstetter, W. Supersolid 
Bose-Fermi mixtures in optical lattices. Phys. Rev. Lett. 
100, 100401 (2008). 

[16] Lewenstein, M., Santos, L., Baranov, M. A. & Fehrmann, 
H. Atomic Bose-Fermi mixtures in an optical lattice. 
Phys. Rev. Lett. 92, 050401 (2004). 

[17] Hubener, A., Snoek, M. & Hofstetter, W. Magnetic 
phases of two-component ultracold bosons in an optical 
lattice. Phys. Rev. B 80, 245109 (2009). 

[18] Capogrosso-Sansone, B., Soyler, §. G., Prokof'ev, N. V. 
& Svistunov, B. V. Critical entropies for magnetic order- 
ing in bosonic mixtures on a lattice. Phys. Rev. A 81, 
053622 (2010). 

[19] Iskin, M. & de Melo, C. A. R. S. Superfluid and insulating 
phases of fermion mixtures in optical lattices. Phys. Rev. 
Lett. 99, 080403 (2007). 

[20] Giinter, K., Stoferle, T., Moritz, H., Kohl, M. & 
Esslinger, T. Bose-Fermi mixtures in a three-dimensional 
optical lattice. Phys. Rev. Lett. 96, 180402 (2006). 

[21] Ospelkaus, S. et al. Localization of bosonic atoms by 
fermionic impurities in a three-dimensional optical lat- 
tice. Phys. Rev. Lett. 96, 180403 (2006). 

[22] Best, T. et al. Role of interactions in ^^Rb-^^K Bose- 
Fermi mixtures in a 3d optical lattice. Phys. Rev. Lett. 
102, 030408 (2009). 

[23] Catani, J., De Sarlo, L., Barontini, G., Minardi, F. & 
Inguscio, M. Degenerate Bose-Bose mixture in a three- 
dimensional optical lattice. Phys. Rev. A 77, 011603 
(2008). 

[24] Weld, D. M. et al. Spin gradient thermometry for ul- 
tracold atoms in optical lattices. Phys. Rev. Lett. 103, 
245301 (2009). 

[25] Gadway, B., Pertot, D., Reimann, R. & Schneble, D. 

Superfluidity of Interacting Bosonic Mixtures in Optical 

Lattices. Phys. Rev. Lett. 105, 21-24 (2010). 
[26] Taglieber, M., Voigt, A.-C., Aoki, T., Hansch, T. W. & 

Dieckmann, K. Quantum degenerate two-species fermi- 

fermi mixture coexisting with a bose-einstein condensate. 

Phys. Rev. Lett. 100, 010401 (2008). 
[27] Taie, S. et al. Realization of a SU(2)xSU(6) system of 

fermions in a cold atomic gas. Phys. Rev. Lett. 105, 

190401 (2010). 

[28] Kitagawa, M. et al. Two-color photoassociation spec- 
troscopy of ytterbium atoms and the precise determi- 
nations of s-wave scattering lengths. Phys. Rev. A 77, 
012719 (2008). 

[29] Batrouni, G. G. et al. Mott domains of bosons confined 
on optical lattices. Phys. Rev. Lett. 89, 117203 (2002). 

[30] Rom, T. et al. State selective production of molecules in 
optical lattices. Phys. Rev. Lett. 93, 073002 (2004). 

[31] Richardson, R. C. The Pomeranchuk effect. Rev. Mod. 
Phys. 69, 683-690 (1997). 

[32] Werner, F., ParcoUet, O., Georges, A. & Hassan, S. R. 
Interaction-induced adiabatic cooling and antiferromag- 
netism of cold fermions in optical lattices. Phys. Rev. 
Lett. 95, 056401 (2005). 

[33] CazaliUa, M. A., Ho, A. F. & Ueda, M. Ultracold gases of 
ytterbium: ferromagnetism and mott states in an SU(6) 
fermi system. New J. Phys. 11, 103033 (2009). 

[34] Hazzard, K. R. A., Gurarie, V., Hermele, M. & Rey, 
A. M. High temperature thermodynamics of fermionic 
alkaline earth atoms in optical lattices. Preprint at 
,http://arxiv.org/abs/1011.0032^ (2010). 



[35] Cramer, M. et al. Do mixtures of bosonic and fermionic 
atoms adiabatically heat up in optical lattices? Phys. 
Rev. Lett. 100, 140409 (2008). 

[36] Fukuhara, T., Sugawa, S., Takasu, Y. & Takahashi, Y. 
All-optical formation of quantum degenerate mixtures. 
Phys. Rev. A 79, 021601 (2009). 

[37] Fukuhara, T., Sugawa, S., Sugimoto, M., Taie, S. & 
Takahashi, Y. Mott insulator of ultracold alkaline-earth- 
metal- like atoms. Phys. Rev. A 79, 041604 (2009). 

[38] Krauth, W., Caffarel, M. & Bouchaud, J.-P. Gutzwiller 
wave function for a model of strongly interacting bosons. 
Phys. Rev. B 45, 3137-3140 (1992). 

[39] Liihmann, D.-S., Bongs, K., Sengstock, K. & 
Pfannkuche, D. Self-trapping of bosons and fermions in 
optical lattices. Phys. Rev. Lett. 101, 050402 (2008). 

[40] Yamashita, M. & Jack, M. W. Mott-insulator shells 
in the three-dimensional bose-hubbard model with har- 
monic confinement. Phys. Rev. A 79, 023609 (2009). 

[41] HackermuUer, L. et al. Anomalous expansion of attrac- 
tively interacting fermionic atoms in an optical lattice. 
Science 327, 1621-1624 (2010). 

[42] Cramer, M. Interaction dependent temperature effects 
in bose-fermi mixtures in optical lattices. Preprint at 
http://arxiv.org/abs/1009.4737 (2010). 



Acknowledgments We acknowledge S. Uetake, T. 
Fukuhara, S. Sugimoto, Y. Takasu and H. Wayama for 
their experimental help and J. Doyle for careful reading 
of the manuscript. This work is supported by the Grant- 
in-Aid for Scientific Research of JSPS (No. 18204035, 
21102005C01 (Quantum Cybernetics)), GCOE Program 
" The Next Generation of Physics, Spun from Universality 
and Emergence" from MEXT of Japan, and World- Lead- 
ing Innovative R&D on Science and Technology (FIRST) . 
S.S and S.T acknowledge support from JSPS. 

Additional Information The authors declare no 
competing financial interests. 



11 



SUPPLEMENTARY INFORMATION 



EFFECT OF THERMAL FLUCTUATIONS IN 
THE REPULSIVELY INTERACTING SYSTEM 

In the repulsively interacting system, we find the 
crossover between the mixed Mott states and the phase- 
separated Mott states. This can be identified by the rapid 
growth of measured bosonic pair occupancy with increase 
in Np as shown in Fig. 3(a). However, from the theo- 
retical curves in Fig. 5(a), one can expect that thermal 
fluctuations easily blur such feature. Our successful ob- 
servation of the crossover is, therefore, due to the cooling 
effect caused by the large entropy carried by fermions as 
mentioned in the main text. Here we discuss in more 
detail how thermal fluctuations affect the crossover be- 
havior of repulsively interacting dual Mott insulators. 

We calculate numerically the temperature dependence 
of the spatial distribution of fermions and bosons. Figure 
[6lja) shows that, in the mixed Mott regime ( Nb = 6x 10'^ 
and Np = 6 x 10^ ), the Mott plateau of total density is 
melted by thermal fluctuations. This leads to the slight 
increase of bosonic pair occupancy with increase in tem- 
perature in this regime as confirmed from Fig. 5(a) at 
Np = 6 X lO'^. On the other hand, in the phase-separated 
regime ( Nb = 6 x lO'^ and Np = 20 x 10'^ ) shown in Fig. 
[S] (b), we find that spatial distributions of bosonic atoms 
become broader and broader with increase in tempera- 
ture, suggesting that the phase-separated state is eas- 
ily destroyed by thermal fluctuations. Consequently, the 
bosonic pair occupancy decreases rapidly with increase 
in temperature as shown in Fig. 5 (a) at Np = 20 x 10"^. 



7|attice=10nK 20 nK 



30 nK 




-40-20 20 40-40-20 20 40-40-20 20 40 

Lattice site 

FIG. 6: Temperature dependence of the spatial dis- 
tribution of fermions and bosons in the repulsively 
interacting system. (a) Mixed Mott regime with A^b = 
6 X 10^ and Nf = 6 x 10"^. (b) Phase-separated regime with 
iVe = 6 X 10^ and iVp = 20 x 10^. Temperature are chosen 



such that Tia 



10, 20, and 30 nK. 



FORMATION OF COMPOSITE PARTICLES IN 
THE ATTRACTIVELY INTERACTING SYSTEM 

In the attractively interacting dual Mott insulators, 
the formation of the composite particles plays a key role 
in the crossover of quantum phases and also thermo- 
dynamics. For further understanding, it is important 
to see how the spatial distributions of such composite 
particles change for different values of Np by carefully 
considering the effect of the trapping potential in the 
experiments. We introduce the local operators A^™™ 
(a =A, B, C and D) that count the numbers of compos- 
ite particles corresponding to type A("BBF"), B("BF"), 
C("BBFF") and D("BFF"] defined in Fig. 4. These 
operators are written as A^^T = |2,1)(2,1|, N^T = 
|1,1)(1,1|, Ng:r = |2,2)(2,2|, and TV^J = |1,2)(1,2| 
(see Method for the notation). The local number oper- 
ators for singly occupied states of bosons and fermions 
are written as TVf'? = |1,0)(1,0| and iVf"? = |0, 1)(0,1|, 
respectively. We calculate the expectation values of these 
operators as a function of Np with fixed A'b = 5 x 10"^ 
and Tini = 40 nK. Figure Ula.) shows the Np depen- 
dence of the total number of the type a composite par- 
ticles, A^^°™ = X]i(^Q°'")' ^i^d singly occupied states, 
iV|" EE (/3 =B or F). In Fig. Wih)-{d), we plot 

the corresponding spatial distributions of singly occupied 
states and composite particles over the sites for several 
Np values. We summarize the formation of composite 
particles in the three different phases discussed in the 
main text; 

• In phase I where Np < 2.5 x 10'^, we see from Fig. 
Ula) that type A("BBF") and type B("BF") com- 
posite particles are created as Np increases. The 
most energetically favorable composite particles of 
type A occupy the sites around the center in the 
trapping potential as shown in the left panel of Fig. 

ETc). 

• In phase II where the number of fermions is inter- 
mediate {Np > 2.5 x 10^), we see in Fig. ^a) and 
(c) that {N^y) and N^^""' grow drastically with 
increase in Np, while {N^"'^) and N'^"^ gradually 
decreases. We also find that the number of the 
bosonic singly occupied state (A^b") decreases sig- 
nificantly in this phase II as shown in Fig. W[.a) and 
(b). In addition. Fig. [TJa) shows that the number 
of composite particles becomes largest in phase II 
around Np ^ 15 x 10^. The similar Np dependence 
is observed for the final temperature Tiattice in Fig. 
5(d) having a broad maximum in phase II. This 
is because the total number of particles effectively 
decreases, and thus the entropy is reduced. 
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FIG. 7: Distributions of the composite particles and singly occupied states, (a) Number of composite particles as 
a function of the number of fermions Nf with A^'b = 5 x 10'^ at Tini — 40 nK. For reference, we also plot the number of singly 
occupied states for bosons (A'^^'"). (b) Spatial distributions of the singly occupied states for bosons (n§" , red dashed line) and 
fermions (iVp", blue line) for several A'^f. (c) Spatial distributions of the composite particles of type A("BBF")((n™™i orange 
line) and type B(("BF")((nB°r, navy dashed line), (d) Spatial distributions of the composite particles of type C("BBFF")((nc,r, 
green dashed line) and type D ( ( "BFF" ) ( (ng"™ , purple hue) . 



• In phase III where Np is large, Fig. Ul^) shows that 
iV^°"^ gradually decreases and both 7V™"^ and iV£°™ 
increase as Np increases. From the right panel of 
Fig. El^d), we find that the more energetically favor- 
able composite particles of type C occupy the sites 
around the center in the trapping potential. In ad- 
dition, the right panel of Fig. ^h) shows that the 
local density of singly occupied state of fermions 
(iVp") grows in the outer region so as to surround 
the type C particles in the central region, showing 
a feature of phase-separation. This originates from 
the effective repulsive interaction between the type 
C particle and the single fermion. Due to the large 
entropy carried by the singly occupied fermions in 
the outer region, the adiabatic heating effects are 
partially compensated in this phase. This results 
in the decrease of the realized Tiatticc in the region 
of A^F > 2.5 X 10" as shown in Fig. 5(d). 



MEASUREMENT OF THE CONTRIBUTIONS 
FROM DIFFERENT COMPOSITE PARTICLES 




5 10 15 20 

A/p [xlO^] 

FIG. 8: Measurement of contribution of composite 
particles of one bosons with one fermions in the at- 
tractively interacting system. Boson-fermion pair occu- 
pancy (red) and boson-fermion pair occupancy after removing 
bosonic double occupancy (green) are measured. From the 
comparison, we can distinguish the contributions from differ- 
ent types of composite particles. The solid lines correspond 
to the numerical simulations for rini=45 nK. 



While we have shown in Fig. 3 the boson-fermion pair 
occupancy measurement in the main text, the measure- 
ments described there cannot distinguish different types 
of composite particles such as type A, B, and C. However, 
it is possible to discriminate the type B ( "BF" ) compos- 
ite particle (single boson and single fermion) from other 
types which has two bosons in addition to a fermion, by 
performing conditional measurements. Figure |5] shows 



two kinds of boson-fermion pair occupancy measure- 
ments at the intermediate atom number of phase II, 
where the type A ("BBF") composite particles are ex- 
pected to be converted into those of type B ("BF"). The 
red curve represents the measurement with the already 
described method. The green one represents the boson- 
fermion pair occupancy measured after removing the sites 
with two bosons by applying the bosonic homonuclear 
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FIG. 9: Number of the composite particles and pair 
occupancies. Various calculated pair occupancies (blue) 
and the associated number of the composite particles (red) 
are plotted as a function of the number of fermions for (a) 
bosonic pair occupancy, (b) fermionic pair occupancy, (c) 
boson-fermion pair occupancy, and (d) boson-fermion pair oc- 
cupancy after bosonic pairs are removed at Tini = 40 uK. 



data of PA measurements. In our experiments, the local 
Fock states consisting of a large number of bosons and 
fermions can be reasonably neglected due to the low fill- 
ing of atoms over the lattice sites. In addition, singly 
occupied states of bosons and fermions are not involved 
in the PA losses. We thus derive approximately the fol- 
lowing useful relations between the pair occupancies that 
are measured by the PA losses and the composite particle 
numbers: 

Db ^ 2Nl°"' + 2iVg°'" , (11) 
Df - 2N^°"' + 2N^°"\ (12) 
Dbf - Nl""^ + TV^"-" + 2Ng"^' + Ar{)°'", (13) 
Dbf-b ~ iV™'" + 7V-"\ (14) 

Here, Dbf-b corresponds to the boson-fermion pair oc- 
cupancy obtained by the conditional PA measurement 
after removing bosonic pairs as explained in the previ- 
ous section. Our numerical results in Fig. ^ support the 
validity of these relations. 



PA laser. The difference in the measured boson-fermion 
pair occupancies corresponds to the contribution from 
the type A ("BBF") and type C ("BBFF") composite 
particles. The measurement clearly indicates that most 
of the composite particles in this regime are of type B 
("BF"), which further supports our explanation. 



EVALUATION OF THE COMPOSITE 
PARTICLES NUMBERS 

It is possible to evaluate the composite particles num- 
bers plotted in Fig. ^a.) directly using the experimental 



PAIR OCCUPANCIES AT VARIOUS LATTICE 
TEMPERATURES 

Here we show in Fig. [10] the numerical calculations of 
bosonic and fermionic pair occupancies as well as boson- 
fermion pair occupancy with various lattice tempera- 
tures. One can also find adiabatic cooling and heating 
effects by comparing these numerical calculations and the 
experimental data shown in Fig. 3(a) and 3(b). 



